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O I In a recent paper, Saxena et al.[l] developed the solutions of three gen- 

^ ■ eralized fractional kinetic equations in terms of the Mittag-Leffier functions. 

■^ . The object of the present paper is to further derive the solution of further 

p ' generalized fractional kinetic equations. The results are obtained in a com- 

pact form in terms of generalized Mittag-Leffier functions. Their relation to 
fundamental laws of physics is briefly discussed. 






Abstract 



1 Introduction 

X. 

H , The fundamental laws of physics are written as equations for the time evo- 

lution of a quantity X (t) , dX (t) / dt = —AX, where this could be Maxwell's 
equations or Schrodinger's equation (if A is limited to linear operators), or it 
could be Newton's law of motion or Einstein's equations for geodesies (if A 
may also be a nonlinear operator [2,3,27]). The mathematical solution (for 
linear operators A) is X{t) = X{0)exp{—At}. 

In thermodynamical or statistical applications one is mostly interested in 
mean values of the quantity < X{t) >. In this case, A is a characteristic time 
scale A~^ = r in the evolution equation for < X{t) >. It then follows that 
< X(t) > decays exponentially toward equilibrium < X{t) >=< (-^(0) > 
exp{— t/r}. 



In 1988, Tsallis [4] generalized the entropic function of Boltzmann-Gibbs 
statistical mechanics, s = — J dxp{x)lnp{x), to nonextensive statistical me- 
chanics with Sg[p] = {1 — J dx\p{x)Y} / {q — 1) that leads to q-exponential 

distributions Pq{x) oc [1 — (1 — q)x'^/kT] . Such a distribution reduces 

to Gaussian distribution for q = 1 and for g = 2 to a Cauchy-Lorentz dis- 
tribution, to name two examples. In an attempt to incorporate Levy dis- 
tribution into statistical mechanics, Tsallis has also shown that the above 
distribution becomes a Levy distribution for q > 5/3. Recently, Tsallis [4] 
used the mathematical simplicity of reaction- type equations, dX/dY = V^, 
to emphasize the natural outcome of the above distribution function Pq{x) 
which corresponds exactly to the solution of the reaction equation of non- 
linear type. The solution has power-law behavior. In the following we show 
that the fractional generalization of the linear reaction-type equation also 
leads to power-law behavior. In both cases, solutions can be expressed in 
terms of generalized Mittag-Leffler functions. 

2 Generalized Mittag-Leffler function 

A generalization of the Mittag-Leffler function [5,6] 

oo n 

EJz) := y — -, (1) 

and its generalized form 

oo n 

EaA^) ■■= T. T^7-—ir^ {a,f3eC,Re{a)>0) (2) 

was introduced by Prabhakar [7] in terms of the series representation 

<.W-|fiit|py (".'5-^^C',fle(„)>0), (3) 

where (7)^ is Pochammer's symbol defined by 

iVn = 1^(^+1). ..(^+„_i)^„gAr 7 7^ (4) 



It is an entire function of order p = [Re{a)]^^ [7]. This function has been 
studied by Wiman [8,9], Agarwal [10], Humbert [11] and Humbert and Agar- 
wal [12] and several others. Some special cases of (3) are given below: 
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eU 



(5) 
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^ , E^,p{z) = Elp{z), (6) 

(m) $(^7;^) = iFMr,z) = T{^)El^{z), (7) 

where $(/5, 7;z) is Kummer's confluent hypergeometric function deflned in 
Erdelyi et al. ([13], p. 248, eq.l]). Mellin-Barnes integral representation for 
the function deflned by (3) follows from the integral 



ElA^) 



2TcujT{'y) Jn 



T{-s)T{j + s){-zy 



Tif3 



sa) 



■ds, 



(8) 



where to = (—1)^/^. The contour fi is a straight line parallel to the imaginary 
axis separating the poles of r(— s) at the points 3 = 1/ (z/ = 0, 1, 2, . . .) from 
those of r(7 + s) at the points s = —7 — u (z/ = 0, 1, 2, . . .). The poles 
of the integrand of (8) are assumed to be simple. (8) can be established by 
calculating the residues at the poles of r(— s) at the points , s = u {u = 
0, 1, 2, . . . , ). It follows from (8) that E^Jz) can be represented in the form 
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I (1-7,1) 
l(0,l),(l-/3,a) 



(i?e(a) >0;a,/3,7GC), (9) 



where Hi2{z) is the H-function. A detailed account of the theory and ap- 
plications of the H-function is available from Mathai and Saxena [14]. This 
function can also be represented by 



ElA^) 



r(7) 



1^1 



(7,1). 

/3,a) ' 



(10) 



where i\E'i(-2) is a special case of Wright's generalized hypergeometric function 

p'^qiz) [15,16] ; also see, Erdelyi et al. ([13], Section 4.1 ), deflned by 



P^9 



\ai,Ai),...,(ap,Ap) 

(bl,Bl),...,{bg,Bg) ' 



UU^{aj + Ajn) z^' 



Snl=ir(6, 



Bju) 
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(11) 



where 1 + Yl'j=i Bj — YTj=i ^j ^ (equality only holds for appropriately 
bounded z). When 7 = 1, (9) and (10) give rise to (12) and (13) given 
below: 
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(13) 



where Re{a) > 0; a, /3, 7 G C. 

If we further take f3 = 1 in (12) and (13) we find that 
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(14) 
(15) 



for Re{a) > 0, z ^ C. The following integral gives the Laplace transform of 
/ e-^'H'^-^El Jar)dt = p-^{l - ap-'')-\ (16) 

JO 
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where Re{p) > \a\ «■=("' , Re{(3) > 0,Re{p) > 0, which can be estabhshed by 
means of the Laplace integral 



e-PHP-^dt 



r(p) 



/o pP 

For 7 = 1, (16) reduces to an elegant formula 



Re{p) > 0, Re{p) > 0. 



/ e-PH'^-^E^ aiandt = p'^^n - ap 
Jo 



a\-l 



(17) 



(18) 
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where Re{j3) > 0,Re{p) > 0,\p\ > |a| «■=(") . In an attempt to investigate 
the functions which when fractionally differentiated (of any order) reappear, 
Hartley and Lorenzo [17] came across a special function of the form 



Fj-a,t] = f^-'Yl 



^0 r(^9 + q) 



Re{q) > 0, 



f^-'E^A-af^). 



(19) 
(20) 



This function has been studied earher by Robotnov [18,19] in connection with 
hereditary integrals for apphcation to continuum mechanics. The Laplace 
transform of this function is given by 

L[Fg{a;t)] = -^,Re{q)>0. (21) 

p'^ — a 

A generalization of the F-function is presented by Lorenzo and Hartley 
[20] by means of the following series representation: 

R.Aa,c,t] = E r[(n+l).-^] '^>->° (22) 

= {t-cy~''-'E,^,_^[a{t-cy],t>c>0. (23) 

The Laplace transform of the R-function is derived by Lorenzo and Hartley 
[20] in the form 

L{K,^{a,c,t)} = ^^, c>0, (24) 

where Rei^u — fj,) > 0, Re{p) > 0. When c = 0, (23) reduces to 

L {R,,^{a, 0,t)} = -^—-, Re{y - /i) > 0, Re[p) > 0. (25) 

Finally we recall the definition of Riemann-Liouville operator of fractional 
integration in the form 

oA-V(t) = TT^ f\t - ur-'f{u)du, (26) 

i [ly) Jo 

with aDif{t) = f(t) [21,22,23]. The standard kinetic equation, when inte- 
grated, yields 

N,it)-No = coD-'N{t), (27) 

where qD^^ is the standard Riemann-Liouville integral operator. Here it 
can be mentioned that in the original paper of Haubold and Mathai [24], 
the number density of the species i,Ni = Ni{t) is a function of time and 
Ni{t = 0) = No is the number density of species i at time t = 0. If we drop 
the index i in (27) and replace c by c*^, then the solution of the generalized 
equation 

Nit)-No = -c''oDrNit), (28) 



is obtained, Haubold and Mathai [24] as 

iVW^iVogyi:!^. (29, 

By virtue of (1), (29) can be written in a compact form as 

Nit) = NoE^i-c'^r). (30) 

In the following, we investigate the solution of two more generalized fractional 
kinetic equations. The results are obtained in a compact form in terms of 
generalized Mittag-LefHer functions and are suitable for computation. A de- 
tailed account of the operators of fractional integration and their applications 
is available from a recent survey paper of Srivastava and Saxena [25]. 

3 Generalized fractional kinetic equations 

Theorem 1. If c > 0, z/ > 0, ;U > 0, then for the solution of the equation 

N{t) - Not'^-'El^l-c'^t^] = -c^oDrNit), (31) 

there holds the formula 

N{t) = Nat^'-^El+\-c''f). (32) 

Proof. By the application of convolution theorem of Laplace transform 
(Erdelyi et al. [26]) we see that (26) can be written as 



p''F{p), (34) 



L[,D;^f{t)-p} = L\^-—\L{m}, (33) 



where F{p) = /o°° e ^^f{u)du,Re{p) > 0. Projecting the equation (31) to 
Laplace transform , we obtain 

[1 + [p/c) ''] 
= Nop-''[l + ip/cyT^^^'^- (35) 
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On using the formula (16), we find that 

L-'ip-'^ {l + ip/cyy^'^'^] = t^-'EI^^-d'f). (36) 

The result (32) now readily follows from (36). 

If we set 7 = 1 then (n)! is cancelled. Then in view of the formula 

PEl,{z) = Ep^,^,{z) + (1 - 7 + P)Ep,,{z), (37) 

which follows from the definition of E2^a{z) given by (3) ,we arrive at 
Corollary 1.1. If c > 0, /i > 0, z/ > 0, then for the solution of 

N{t) - N^t^-^E,,^[-cH-] = -d' oD^Nit), (38) 

there holds the formula 

N{t) = -^[E,,,^,{-c''tn + (1 - /i + u)E,,,{-c''n]. (39) 



^^,7^ = 77J^[EprMz) + {3P + 3-2^)Ep,^^,{z) 



When 7 = 2, then by virtue of the following identity 

1 

2P' 

+ {2/32 + 72 + 3/?-27-3/37 + l}^A7(^)], (40) 

which follows as a consequence of the definition (3) , we obtain 
Corollary 1.2. If c > 0, z/ > 0, /i > 0, then for the solution of 

N{t) - Nor-'ElJ-c-'n = -c^^D^Nit), (41) 

there holds the relation 

N{t) = Not^-'El^i-c'^n, (42) 

= -^ [E,,,-2{-c''n + {3(z/ + 1) - 2fi} E.^.^.i-c'tn 
+ {2(z/2 + ^2) + 3z/-2/i-3z//i+l}E,,^(-c"0]. (43) 

Next, if we set 7 = 0, then by virtue of the identity E^^^i^z) = pT-r, we arrive 
at another result given by Saxena et al. [1]. 



Theorem 2. If c > 0, 6 > 0, Re{p) > 0, z/ > /i + 1, then for the solution of 

N{t) - NoR^^^i-c", b, t) = -c" oD^'^Nit), (44) 

there holds the formula 



iVf 



N{t) = -^{t - 6)^-^-Mi?.,.-^-i(-c^(t - by) + (/i + i)i?.,,_^(-c'^(t - by 



Proof. Taking Laplace transform of both sides of (45), it gives 
N{t) = L{N{ty,p} = NoL-' 



(45) 



^-bp fi-u 



{1 + (c^/p-)}' 



N^L-^ 



-bp fi—u 



Y"' E 
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iVoE 
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[^)n[ C) j--lr -bp u-u-nu-i 



^0 {ny.T{nu + u - fi) 
-^ ^ ^0 r(nz/ + z/-/i) 
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iVo 
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iVo, 



it-b) 



u—fi—l 






_^(t _ 6)-^-i[i?,,,_^_i {-c-(t - by] + (^ + i)E,,,_^ {-c-(t - bf]], 

which is same as (45). This completes the proof of theorem 2. If we set 
/i = 0, theorem 2 reduces to 

Corollary 2.1. If c > 0, 6 > 0, z/ > 1, then for the solution of 

N{t) - iVoi^.,o(-c^ &, t) = -c" oD^'Nit), (46) 

there holds the formula 



N{t) 



Nn 



V 



(t - br-'[E,,,_,{-c''{t - bY) + E,4-c''{t - byy. 



(47) 



For 6 = 0, theorem 2 yields 



Corollary 2.2. If c > 0, z/ > /i + 1, then for the solution of 

iV(t) -iVoi?,,^(-c^O,t) = -d" oD^^Nit), (48) 

there holds the formula 

N{t) = ^t''-^^-\E,,u^^_,{-c''n + (^ + l)E,,,_^(-c-r)]. (49) 

If we further take /i = then the above corollary reduces to the following 
result: 

If c > 0, z/ > 1, then the solution of 

N{t) - A^oF,[-c^ t] = -c'D^'Nit), (50) 

is given by 

Nit) = ^-—[E,^,_,(-c''tn + E^A-c'^n]- (51) 
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